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We present a model describing evolution of the small-scale Navier-Stokes turbulence due 
to its stochastic dis tortions by much la rger turbulent scales. This study is motivated by 
numerical findings I Laval et all l|200ll) ) that such interactions of separated scales play 



important role in turbulence intermittency. We introduce description of turbulence in 
terms of the moments of the fc-spa ce quantities using a m ethod previously developed 
for the kinematic dynamo problem l|Nazarenko et jjD ihuO.Sh ). Working with the fc-space 
moments allows to introduce new useful measures of intermittency such as the mean po- 
larization and the spectral flatness. Our study of the 2D turbulence shows that the energy 
cascade is scale invariant and Gaussian whereas the enstrophy cascade is intermittent. 
In 3D, we show that the statistics of turbulence wavepackets deviates from gaussianity 
toward dominance of the plane polarizations. Such turbulence is formed by ellipsoids in 
the fc-space centered at its origin and having one large, one neutral and one small axes 
with the velocity field pointing parallel to the smallest axis. 



1. Introduction 

Finding a good turbulence model is a long standing problem. To be useful in appli- 
cations, the model has to be sufficiently simple and yet capable of capturing the basic 
physical processes such as the energy cascade and the intermittent bursts. The cascades 
appear to be a more robust property reasonably well descri bed by classical tu rbulence 
closures such as the direct interac t ion ap proximation (DIA) l|Kraichnanl l|l96lf) ) and its 
derivatives (e.g. EDQNM lOrszaj (^)66)). The turbulence intermittency appears to be 
a more subtle process which depends on the detailed features of the dynamical fluid 
structures. Of particular importance is the question whether the intermittent bursts are 
caused by finite-time vorticity "blow-ups" (believed to become real singularities in the 
limit of zero viscosity) or a "slower" exponential vortex stretching by a large-scale strain 
collectively produced by the surrounding vortex tubes. Note that the first process is lo- 
cal in the scale space, - it is usually viewed as two or more vortex tubes of similar and 
implosively decreasing radius. On the other hand, the second process involves interaction 
of significantly sepa rated scales, - a thi n vortex tube and a large-scale strain. Recent nu- 
merical simulations l|Laval et all l|200ll) ) indicate that it is the nonlocal scale interactions 
that are responsible for the deviations of the structure functions from their Kolmogorov 
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self-similar value whereas the net effect of the local interactions is to reduce these devi- 
ations. These conclusions lead to a model of turbulence in which a model closure (e.g. 
DIA) is used for the local interactions whereas the nonlocal interactions are described by 
a wavepacket (WKB) formalism which exploits the scale separation. The later describes a 
linear process of distortion of small-scale turbulence by a strain produced by large scales. 
Such a linear distortion is a familiar process in engineering applications, for example, 
when a turbulent fluid flows through a pipe with a sudden c hange in diameter. It is 
descri bed by the rapid distortion theory (RDT) introduced bv iBatchelor k. Proudmanl 
The model considered in this paper is different from the classical RDT in that 
the distorting strain is stochastic and, therefore, we will call it the stochastic distortion 
theory (SDT). We will study a simplest version of SDT in which the large-scale strain is 
modeled by a Gaussian whit e (in time) noise i n the spirit of the Kraichnan model used 
for turbulent passive scalars llKraichnanI dl974ln and of the Kazantscv-Kraichnan model 
from the turbulent dynamo theory IjKazentsevI l|l968|) : Fkraichnan fc Nagaraianl l)l967(l ). 

Because most of studies so far have focused on the theories with local scale interactions, 
we will devote our attention mainly to the description of the nonlocal interactions. The 
local interactions are unimportant at small scales in 2D, but they should be taken into 
account when 3D turbulence is considered. Similarly to RDT, the SDT model deals 
with fc-space quantities due to much greater simplicity of the pressure term in the fc- 
space. We will study statistical moments of the fc-space quantities of all orders and 
not only the second order correlators as it is customary for RDT. The higher fc-space 
correlators carry an information about the turbulence statistics and intcrmittency which 
is not always available from the two-point coordinate space correlators, the structure 
functions, which are popular objects in the turbulence theory. Indeed, intermittency in 
some systems can be dominated by singular fc-space structures which are not singular in 
the i-space (e.g. periodic fields). These structures leave their signature on the scalings of 
the fc-space moments but not on the x-space structure functions. The system considered 
in the present paper is of this type, and anoth er example of thi s kind is the magnetic 
turbulence in the kinematic dynamo problem l|Nazarenko et a,l,\ ty2QQ?t )). In fact, SDT 
bears a lot of similariti es to the turbulen t dyna mo problem and in this paper we will use 
a method developed in lNazarenko et al\ l|2003f) for its derivation. We will also see that, 
like in the dynamo problem, the fourth order fc-space moments allow us to introduce the 
measures of the mean polarization and of the spectral flatness, - the quantities of special 
importance for characterization of the small-scale turbulence. 

2. Stochastic distortion of turbulence 

Let us consider a velocity field in three-dimensional space that consists of a component 
U with large characteristic scale L and a component u with small characteristic scale Z, 
L ^> I. In this case, Navier-Stokes equation is 



d t U + <9 t M+(U-V)U+(U- V)u + (u- V)U+(u- V)u = -Vp + vV 2 V + vV 2 u. (2.1) 

Let u s define the Gabor tran sform (GT) (see iNazarenko fc Laval l)200f)lk iNazarenkol 
(1999'l: lNazarenko et oil j200C|U 



where 1 » f ' > e and f(x) is a function which decreases rapidly at infinity, e.g. exp(— x 2 ). 
Averaging (•) is performed over the statistics of a random force which will be introduced 
below. 




(2.2) 
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One can think of the GT as a local Fourier transform taken in a box centered at x 
and having a size which is intermediate between L and I. The GT commutes with the 
time and space derivatives, dt and V. Commutativity with dt is obvious. Note that the 
GT commutes with V only for distances from the boundaries which are larger than the 
support of function /. The inverse GT is simply an integration over all wavenumbers, 
e.g. 

If dk 
u ( x > t ) = J^Ju( x > k > t )j2^3- ( 2 - 3 ) 

Here, we will study only the nonlocal interaction of small and large scales and therefore 
we neglect the nonlinear term (u-V)u which corresponds to local interactions among the 
small scales. Let us apply the GT to the above equation with k ~ 2n/l ~ 1 ^> 2tt/ L ~ e 
and only retain terms up to first power in e and e* (we chose e* such that £*>e> ( e *) 2 )- 
All large-scale terms (the first and the third ones on the LHS and the third one on the 
RHS) give no contribution because their GT is e xponentially sma l l. Equ ation for the GT 
of u under such assumptions where obtained in lNazarenko et all l|2000(l : it is 

2k 

D t u+(u-V)V = -rj u-V(\J -k)-vk 2 u, (2.4) 

where 

D t = d t + x- V + k V k , 
± = U, (2.5) 
fe = -V(fe-U), (2.6) 

Equation (|2.4|l provides an RDT description of turbulence generalized to the case when 
both the mean strain and the turbulence are inhomogeneous. This equation has the form 
of a WKB-type transport equation with characteristics given by 1)2. 5fl and 1)2. 6f> . Consider 
this equation for a fluid path determined by x(t) = U, so that u(k, x, t) — > u(k, x(t), t) f 

2 

dtu m = o~ijki djU m o~ m iUi + -j^k m (o~ijkiUj) vk u m , (2-7) 

where o~ij = V jU% is the strain matrix and operators Vj and di mean derivatives with 
respect to Xi and fc, correspondingly [i — 1, ...,£>). Note that strain (taken along a 
fluid path) enters this equation as a given function of time. Equation (|2.7|l is applicable 
to arbitrary slowly varying in space large-scale flow. We formulate the SDT model as 
equation l|2.7(l complemented by a prescribed statistics of the large-scale flow. One can 
use, for example, a numerically computed large-scale strain and use it as an input into the 
equation (|2.4(l should later be integrated numerically. In this paper, however, we would 
like to derive a reduced model via the statistical averaging which is possible by assuming 
a sufficiently simple statistics of the large-scale strain. Experiments and numerical data 
indicate that Navier-Stokes turbulence is Gaussian at large scales and we will use this 
property in our model . We will further assume that the strain is white in time with 
iNazarenko et al\ l|2003|) 

try = - ^%) (2.8) 

where Aij is a matrix the elements of which are statistically independent and white in 



f Hereafter, we drop hats on u because only Gabor components will be considered. Also, we 
will not mention explicitly dependence on the fluid path and simply write u = u(k, t). 
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time, 

(Aij(t)Aki(0)) — 5ijSki5(t). (2.9) 
This choice of strain ensure the incompressibility and statistical isotropy. In this case 

K (iH,(0)) = fl(S ik 6ji - ^SijSkOSit). (2.10) 

Note that this is not the only way to satisfy the incompressi bility and the i s otrop y 
and there are infinitel y many ways to choose such s tatistics; e.g. IChertkov et oil (^999), 
iFalkovich eTaW l|200lh . iBalkovskv fc Fouxonl l|l999h choose 

(ffii(t)o*i(0)) = n((d + ySikdfl - SaSjk - 5 t3 5 kl )S(t). (2.11) 

However, any such choice will lead to (up to a time-scale constant) the same final equa- 
tions. The Gaussian white in time strain has a lso be e n used in the MHD dynamo the - 
ory (Kazantsev-Kraichnan model: iKazentsevI Jl96ST> . iK raichn an k. Naearaianl l|l967f) ) 
and in the theory of turbulent passive scalar l|KraichnaiT " i}l974|k see also review of 
Falkovich et al It is a natural starting point because of its simplicity. 

Note that for realistic modeling of small-scale turbulence one has to describe a matching 
to the large-scale range via a low-fc forcing or a boundary condition. As we will see 
later, some properties of the small-scale turbulence turn out to be independent of these 
effects, e.g. scalings in 2D case, whereas the 3D case is more sensitive to the boundary 
conditions. Detailed modeling of the \ow-k forcing/boundary conditions is beyond the 
scope of this paper. Below, we will simply consider forcing-free evolution of a finite- 
support initial condition (decaying turbulence) . We will also consider finite-flux solutions 
in 2D corresponding to the turbulent cascades in forced turbulence. 



3. Generating Function 

Let us consider the following set of 1-point correlators of Gabor velocities, 

^ = (|u(k)| (2 "- 4s) |u(k)T s > (3-1) 

with n = 1, 2, 3, ... and s = 0, 1, 2, 3, .... Such correlators where shown in lNazarenko et all 
( 2003) to be a fundamental set in case of homogeneous isotropic turbulence from which 
one can express any of two-point correlators of the following kind 

(ui 1 (ki)ui 2 {ki)..M ln (k 1 )u jl (k2)u j2 (k2)...u Jm (k 2 )), (3.2) 

where ii,i2, ■■,in and ji,j2, ■-,jm take values 1, 2 or 3 indexing the components in 3D 
space, and n and m are some arbitrary natural numbers. Note that homogeneity and 
isotropy in SDT follow from the coordinate independence and isotropy of the strain and 
it has to be understood only in a local sense, near the considered fluid particle of the 
large-scale flow. 

Now we define a generating function, 

Z(X, a, (3, k) = (eAluO^+auO^+zJudc)^ (3 3) 

where overline denotes the complex conjugation. This function allows one to obtain any 
of the fundamental 1-point correlators H3.1(l via differentiation with respect to A, a and 

(3.4) 

To derive an evolution equation for Z we will follow the technique developed in 



d 



(2ra-4s) 
A 



d» a d s Z 
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iNazarenko et all (|2003h for the turbulent dynamo problem. Let us time differentiate 
the expression for Z (|3.3|1 and use the dynamical equation (|2.7(l : we have 

Z = k l dj{o- l jE) - \{a m i(u m ui +uiu m )E) - 2a{a m iu m uiE) - 2f3{a m iu m uiE) 

-2wfc 2 ((A|u(k)| 2 + au(k) 2 + (M(k) 2 )E) (3.5) 

where 

E = e A|u(k)| 2 +Qu(k) 2 +/3u(k) 2 ^ /ggS 

To find the correlators on the RHS of l|3.5[l . we use Gaussianity of cry and perform a 
Gaussian integration by parts. Then, we use whiteness of the strain field to find the 
response function (functional derivative of ui with respect to cry). Finally, we use the 
isotropy of the strain so that the final equation involves only k — |k| and no angular 
coordinates of the wave vector. Leaving the derivation for the Appendix, we write here 
only the final result, 

z = n[(i- \)k 2 z kk + hiv + d 2 - i)k,z k + 2(i ~ ^- + d)vz - ^v 2 z 

+2(A 2 + 4a[3)Z xx + 2\ 2 Z al3 + 8XaZ aX + 8Xf3Z fjx + Aa 2 Z aa + 4/3%^] 
-2vk 2 VZ, (3.7) 

where the k,a,(3 and A subscripts in Z denote differentiation with respect k,a,(3 and A 
correspondingly and 

V = Xd x +ad a +/3d/3. (3.8) 

The number of independent variables in this equation can be reduced by one taking into 
account that due to turbu l ence h omogeneity Z depends on a and j3 only via combination 
j] = a/3 l|Nazarenko et al\ ll2003[) 1. We have 

Z = 0[(1 - ^)k 2 Z kk + i(4£> + d 2 - l)kZ k + 2(1 - ^ + d)VZ - ^D 2 Z 

+2(A 2 + 4?7)Z AA + 2\ 2 {Z n + r)Z nn ) + W\r)Z nX + &v 2 Z vv ] - 2vk 2 VZ, (3.9) 

where 

V = Xd x + 2r)d v . (3.10) 

Equation (|3.9|l is the main equation of SDT. The RHS of this equation describes inter- 
actions of the separated scales only. In practical applications or numerical modeling one 
has to add to it a suitable model for the local scale interactions. We leave this task for 
future and concentrate below on studying the effect of the nonlocal interactions only. 



4. 2D turbulence 

Let us first of all consider the 2D case. The large time dynamics of the small-scale 
turbulence is known to be dominated in 2D by the nonlocal interactions due to generation 
of intense large-scale vortices and, therefore, the SDT model l|3.9|l is relevant even without 
including a model for the local interactions. Note that equations for the 2D turbulence 
in which only nonlocal interactions are left are formally identical to the passive scalar 
equations. The energy spectra of the nonlocal 2D turbulence an d the passive scalars 
in the Batchelor regime were studied in INazarenko fc Laval (|2000t) without making any 
assumptions on the strain statistics. Here we will study the higher fc-space correlators. 

The 2D case is simpler than the 3D one in that all correlators 'J'™ with s > are not 
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independent and can be expressed in terms of 'Fq which we will call the energy series, 

= E n (k, t) = (|u(k)| 2 ") = [d^Z] x=n=0 . (4.1) 

The equations for correlators E n are to be obtained by differentiating l|3.9[l n times with 
respect to A and taking it at 77 = A = which gives 

E n = - [k 2 {E n ) kk + (3 + An)k{E n ) k + 4n(l + n)E n ] - 2vnk 2 E n . (4.2) 

Of special interest are the cascade-type solutions realized when turbulence is forced. To 
obtain these solutions we hrst re-write 1)4. 2f) as a continuity equation in fc-space; this can 
be done in two different ways, 

dtik^En) = ^[k- 1 (k 2n+2 E n ) k ] k - 2vnk 2n+1 E n . (4.3) 

and 

d t {k 2n+1 E n ) = ^[k 3 (k 2n E n ) k ] k - 2vnk 2n+3 E n . (4.4) 

In the absence of dissipation, v = 0, equations 1)4.4(1 and (|4.3(l describe conservation of the 
quantities which have spectral densities £ n = k 2n ~ 1 E n and T n = k 2n+1 E n respectively. 
For n — 1 these quantities are just the energy and the enstrophy f . We will call invariants 
£ n and T n the energy and the enstrophy series correspondingly. The steady state solutions 
in the range where viscosity is negligible are 

E n = C[ n) k~ 2n (4.5) 

and 

E n = C ( 2 n) k- 2n - 2 , (4.6) 

where c[ n ^ and are arbitrary po sitive constants. For n — 1 these solutions where 
obtained in lNazarenko fc Lavall l)2000(k solution l|4.5|l corresponds to a constant energy 
flux (and equipartition of the enstrophy) whereas (|4.6() corresponds to an enstrophy 
cascade (and equipartition of the energy). Because the equations for E n are linear, any 
linear combination 

E n = C[ n) k- 2n + C ( 2 n) fc" 2 "" 2 (4.7) 

is also a stationary solution (in fact the general one) . Note that the flux of invariant £ n is 
given by —^k~ 1 {k 2n+2 E n ) k and it is always negative on solutions l|4.7(l whereas the flux 
of T n , which is —^k 3 (k 2n E n ) k , is always positive in the steady state. Thus, in forced 
turbulence solutions (|4.5|l will form on the low-fc side of the forcing scale and solutions 
(|4.6|) on the high-fc side of it, which agrees with the general observation that the energy 
cascade is inverse and the enstrophy cascade is a direct one. 
Let us introduce a spectral flatness, 

F n = E n /E%. (4.8) 

For Gaussian fields, F n would be independent of fc and, therefore, fc-dependence of F n 
bears an information about the scale invariance and presence of turbulence intermittency. 
In particular, for the solutions (|4.5I) we have F n — const indicating that the inverse 

f Recall that we consider nonlocal turbulence and, although invariance of enstrophy is easily 
seen from the vorticity conservation along the fluid paths, the energy invariance is not obvious 
because there can be non-local energy exchanges between turbulence and the large-scale flow. 
The c onservation of energy was proved for initially isotropic turbulence in Nazarenk o fc Lavall 
(2000) 
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cascades are not intermittent: turbulence produced by a Gaussian forcing at some scale 
kf will remain Gaussian at k < kf. On solutions (|4.6(l we have F n ~ fc 2n ~ 2 which indicates 
broken scale invariance and growing deviation from Gaussianity at small scales. Such a 
small-scale intermittency in the direct cascades is due to nearly singular fc-structures 
having the shape of strongly elongated ellipses in 2D fc-space which are centered at the 
origin. Each strain realization will produce just one of these (randomly oriented) ellipses 
out of an initial circular (isotropic) distribution. 



5. Nonlocal 3D turbulence 

3D Navier-Stokes turbulence is never likely to be nonlocal and for its realistic de- 
scription one should add a model of the local interactions into the RHS of equation 
l|3.9|> . However, it is still interesting to solve equation l|3.9[l as is in order to study the 
effect of pure nonlocal interactions. We will see below that such a study will reveal some 
interesting physics. 

In the 3D case, equation 1)3.9(1 is 

Z=— [k 2 Z kk + AkZ k + (8 + 2kd k )VZ + {\ 2 Z XX + 8r 1 2 Z m ) 

+ (3A 2 - 4jj)(Z„ + T]Z m ) + l6\vZ Xr , + l2r,Z xx ] - 2vk 2 VZ, (5.1) 

The standard procedure to obtain equations for the correlators is to differentiate 
l|5.1|) with respect to A and rj the required number of times and then to put A = rj = 0. 
In the 3D case, all the correlators \&™ are independent and at each order 2n we have a 
system of coupled equations rather than a single equation to solve as it was the case in 
2D. However, a decoupling arises asymptotically at large times as we will see below. We 
will start by considering the second and the fourth order correlators (n — 1 and 2). 

5.1. Energy spectrum 

E(k, t) ee E 1 (k, t) = (|u(k)| 2 ) = [d x Z] x=v=Q . (5.2) 
Differentiating (|5.1() with respect to A and taking the result at A = r\ = we have 

20 

E = — (k 2 E kk + 6kE k + 8E) - 2vk 2 E. (5.3) 

This equation is similar to the Kazantsev equation ijKazentsevI (|l968^) ^1 describing evolu- 
tion of the magnetic energy spectrum in the kinematic dynamo theory. Similarly to the 
dynamo theory, the total energy grows exponentially and therefore, unlike the 2D case, 
no stationary cascade states are possible. To have a steady state, one has to add a model 
of local interactions to SDT which will be done in future publications. However, we will 
now study equation (|5.3I) to examine consequenc es of interactions of separat ed scales. 

Recently, Schekochihin, Boldyrev and Kulsrud l|Schekochihin et~al\ l|2002tJl ) presented 
the solution of the Kazantsev equation obtained by the Kontorovich-Lebedev transform 
and we will use their results to solve 1(5. 2J) . By substitution 



E = e int ' 6 k- 5 ' 2 (j ) {k/k cU t) k d = x /n/3^, (5.4) 
one can reduce l|5.3|l to 

3 • 

■^(p{p,t)^p 2 (j)pp+p<j)p-p 2 <j). (5.5) 

The RHS of this equation is just the modified Bessel operato r and by using the Kontoro vich- 
Lebedev transform one immediately gets for t ^> 1/fi: IjSchekochihin et all ((2002^ . 
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iNazarenko et all ^200^) ) 

0(p, t) = const / ds s sinh(7rs) K is (p)K is (q) e^* 2 ™/ 3 , ( 5-6 ) 
Jo 

where Ki S is the MacDonald function of an imaginary order and the constant is fixed by 
the initial condition. At scales much greater than the dissipative one, p <C 1, viscosity is 
not important for t <C (lng) 2 ( q <C 1 is the mean wavenumber of the initial c ondit ion) 
and the solution (|5.6|) becomes l|Schekochihin et all l)2002aj) ; INazarenko et all l|2003fl 'l 

, ,„ 3(lnfc/g) 2 

4> = const t~ 1/2 e Brii — , (5.7) 

This solution describes a spectrum with an expanding fc~ 5 / 2 scaling range (which means 
fc~ 1//2 for the one dimensional energy spectrum). At t ~ (hit?) 2 the front of this scaling 
range reaches the dissipative scales and for t ^> (mg) 2 and f 15.60 gives 

ct,= C ^K (p), (5.8) 

Function Kq(p) decays exponentially at large p which corresponds to a viscous cut-off of 
the spectrum. For p <C 1, Kq(p) w — lnp, which means that at large time the scales far 
larger than the dissipative one are affected by viscosity via a logarithmic correction, 

E(k) = ^ e ™'/ 6 fc- 5 / 2 ln(Vfc)- (5.9) 



6. 4th-order correlators, turbulence polarization and flatness 

There are two independent 4th order correlators 

S(k,t) = ^ 2) = (|u(k)| 4 > = \Z XX ] X=V=Q and T{k.t) = ^ = (|u 2 (k)| 2 ) = [Z,] x= ^ 

(6.1) 

Differentiating l|5.1|l twice with respect to A and taking the result at A = r\ = we have 

20 

S = —(k 2 S kk + 8kS k + 185 + 6T) - 4i/fc 2 5. (6.2) 
Now, differentiating l|3.9|l with respect to r\ and taking the result at A = r\ = we get 

T = — {k 2 T kk + 8kT k + 12T + 125) - \vk 2 T. (6.3) 
o 

By subtracting (|6.3|l from 16.211 we get a closed equation for W = 5 — T, 

on 

W= — (k 2 W kk + 8kW k + 6W) - 4vk 2 W. (6.4) 
The physical meaning of W becomes clear if we re- write it as INazarenko et all |2003) 

3 3 

W = A^m^i)} 2 ) = 4]T(K| 2 M 2 sm 2 (0 J - ^)> > 0, (6.5) 

where symbol 9 denotes the imaginary part and <f>j and 4>i are the phases of compo- 
nents Uj and ui respectively. Thus, we see that W contains information not only about 
the amplitudes but also about the phases of the Fourier modes. In particular, 14^ = 
corresponds to the case where all Fourier components of the magnetic field have plane 
polarization. If W ^ then other polarizations (circular, ecliptic) are present. This is the 
case for example for Gaussian fields when W = E 2 /2 > 0. On the other hand, smallness 
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of the phase differences can be overpowered in W by large amplitudes. Therefore, a better 
measure of the mean polarization would be a normalized W , e.g., 

P = W/S. (6.6) 

Defined this way the mean turbulence polarization is an example of an important physical 
quantity which can be obtained from the one-point Fourier correlators and is unavailable 
from the coordinate space (one-point or two-point) correlators. 

Equation l|6.4|l can be solved similarly to the energy spectrum equation i|5.3|) . namely 
by transforming it into equation l|5.5|l by a substitution similar to (|5.4|) and then using 
the solution (|5.6|l . In the inviscid regime ( p <C 1, t <C (\nq) 2 ) we have 

^ = V^ ^ 1/2 e- 25f2t/6 fc- 7/2 e" H1 ^, (6.7) 

where Wo is a constant which can be found from the initial conditions. We see that 
W develops a fc~ 7 / 2 scaling range which is cut off at low and high k by exponentially 
propagating fronts. Within this scaling range, W decays exponentially in time. 

Given W, one can find S by using substitution S = V + W/3 which leads to a closed 
equation for V which can be solved similarly to E and W. This gives for the inviscid 
regime 

S = r 1 ' 2 k- 7 / 2 e- 3J ^- (Vo e 47nt / 6 + ±W e ~ 25nt ^ , (6.8) 

where Vq is another constant which can be found from the initial conditions. For i > 1, 
the second term in the parenthesis should be neglected. Then, we have the following 
solution for the mean turbulence polarization, 

P = W/S= — E~ UQt . (6.9) 
Vo 

As we see, in the inviscid regime the mean polarization tends to an independent of k 
value which exponentially decays in time. This means that all turbulence wavepackets 
eventually become plane polarized. Recall that such turbulence is very far from being 
Gaussian for which the mean polarization remains finite (elliptic and circular polarized 
modes are present). 

In the diffusive regime, t 3> (lnq) 2 , we have 

W(k) = W i- 1 /2 e -250t/6 k -7/2 fc^/y^.) ( 6 . 10 ) 

and 

S(k) = Vo r^e 471 */ 6 k- 7 ' 2 \n(k d /V2k). (6.11) 

Thus, P is still give by the same formula l|6.9fl indicating that the mean polarization 
continues to further decrease in time with the same exponential rate. Thus, by the time 
the diffusive regime is achieved W can be essentially put equal to zero. 

The fact that the polarization becomes plane has quite simple physical explanation. 
Indeed, a vorticity wavepacket of arbitrary polarization will be strongly distorted by the 
stretching which mostly occurs along the dominant eigenvector of the lagrangian defor- 
mation matrix (corresponding to the greatest Lyapunov exponent). Such a stretching 
make any initial "spiral" flat for large time with the dominant field component lying in 
the plane passing through the stretching and the wavevector directions (and, of course, 
u(k) is perpendicular to k). 

Another important measure of turbulence intermittency available from the A:-space 
moments is the spectral flatness which can be defined as F — S/E 2 . For large time in 
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the inviscid regime 

i^l/2 e -lim/2 fc 3/2_ (gl2) 

We see that the flatness growths both in time and in k which indicates presence of the 
small-scale intcrmittency. Such an intcrmittency can be attributed to the presence of 
coherent structures in fc-space. 

One can also find solution for the fourth order correlator S in the dissipative regime. 
This will be done in the next section together with correlators of all higher orders. 



7. Large-time behavior of higher correlators 

The observation in the end of the previous section (that there is a dominant field 
component) allows us to predict that for large time |u| 4 « |u 2 | 2 in each realization so 
that Z\\ ~ Z a @. Therefore, property Z\\ = Z a p, if valid initially, should be preserved 
by the equation for Z. Indeed, by using the equation (|5.1|l combination w = Z\\ — Z a p = 
Z\\ — Z n — rjZ^ satisfies a closed homogeneous equation. This means that if w = at 
t = then it will remain identically zero for any time. Thus, we can consider a class of 
solutions of (|3.9[1 (corresponding to large-time asymptotics of the general solution) such 
that Z\\ = Z a p. Assuming this equality in 13.9fl and putting r\ — we have, 

Z=— [k 2 Z kk + 4kZ k + (8 + 2kd k )XZ x + 4A 2 Z AA ] - 2^ 2 AZ A , (7.1) 

This gives the following equations for the correlators E n = (|u(k)| 2 ™) is the correlator of 
order 2n, 

E n = — [k 2 {E n ) kk + (2n + A)k{E n ) k + An(n + l)E n ] - 2isk 2 nE n . (7.2) 

Note that for n = 1 this coincides with for the energy spectrum equation (|5.3(l which we 
have already solved. Moreover, by substitution 

E n = e (3« 2 +«-9/4)2Qt/3 k (-n-3/2) 0( fc ^/^ t) (7.3) 

one can reduce (|7.2() to an independent of n equation (|5.5() for function 4> the solution of 
which we already know. For the inviscid regime (I/O <C t <Si (Inq) 2 ) we have 

E n - COnst(Tl) r 1 ' 2 e (3n 2 +™-9/4) 2 O t /3-Hl^ fc (-n-3/3) j (7.4) 

and in the diffusive regime (t 3> (lng) 2 ) 

En = COnS t(") e (3n' + n-9/4)2n t /3 fc («-3/2) K ^ k/kd y (7.5) 

This expression agrees with the results obtained in the previous sections for n = 1 and 
n = 2. We see that the main effect of the dissipation is the prefactor change t~ x l 2 — * 
i~ 3 / 2 and the K§(kjkd) form-factor which corresponds to a log-correction at k <C 1 and 
an exponen tial cut-off at k ~ fa. Similar results were obtained for the magnetic field 
moments in lNazarenko et all l|2003|) . In that case it is the exponential cut-off that cause s 
the change of the exponent ial growth in the mean magnetic energy ijKazentsevI l(l9681. 
iKulsrud fc Anderson! |l9 92'0 and in the higher x-space moments of the magnetic field 
l|Chertkov et all l|l99fl) ~ 

The scalings in l|7.4|l and (|7.5() with respect to the order n contain important informa- 
tion about the small-scale turbulence. In particular the exponential growth in time with 
exponent ~ n 2 indicates that the turbulence statistics is log-normal. An equivalent result 
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for the mag netic fields in the kinematic dynamo problem was obtained in lChertkov et all 
l|l999h and ISchekochihin et all [|2002^ . The log-normality arises because the strain is 
a multiplicative noise for the velocity field which becomes nearly one-dimensional and 
because the time integrated strain tends to become a Gaussia n process. Fo r mally , this 
resu lt can be also obtained using the random matrix theory of iFursteneerd l)l963[) (see 
also iFalk ovich e t all (l2001^l > ). and a more detailed physical explanation can be found in 
iNazarenko et all ll2003j) ~ 

The fc-dependence of the Fourier correlators is also very important because it gives 
an information about the dominant structures in the wavenumber space. Suppose that 
initially the turbulence is isotropic and concentrated in a ball centered at the origin in 
the wavenumber space. For each realization such a ball will stretch into an ellipsoid with 
one large, one short and one neutral dimensions. One can visualize this ellipsoid as an 
elongated flat cactus leaf with thorns showing the velocity field direction. Note that in 
this picture one component of the velocity field (transverse to the cactus leaf) is dominant 
which is captured by the fact that the polarization W introduced in this paper tends to 
zero at large time. Another consequence of this picture is that the wavenumber space 
will be covered by the ellipsoids more sparsely at large k which implies large intermittent 
fluctuations of the velocity field in the /c-space. These fluctuations could be quantified 
by the flatness F which was shown in (|6.12|) to grow as A; 3 / 2 , a clear indication of the 
small-scale intermittency. 

8. Sensitivity of SDT to the strain statistics 

In the previous sections, the SDT model was formulated and studied assuming that 
the large-scale strain is a Gaussian white noise process. This assumption, similar to the 
Kazantsev-Kraichnan dynamo model, allowed us to obtain some important analytical 
solutions for the energy spectrum, the polarization, the flatness and higher-order corre- 
lators which capture the turbulence intermittency. In real experimental and in numerical 
turbulence the strain is generally far from being a Gaussian white noise. Thus, it is inter- 
esting to study sensitivity of our SDT model and its predictions to the strain statistics. 

In this section, we will numerically simulate l|2.7|) for two different types of strain. 
First, we consider a synthetic Gaussian strain field with a finite correlation time r which 
is algorithmically generated as 



where Aij is the same matrix as in (|2.9|) and dt is the time step. The r.m.s. of the different 
strain components in this numerical experiment ranged from 2.9 to 3.6 and the correlation 
time was 0.02. Thus, the correlation time was about 16 times less than the characteristic 
strain distortion time. We considered 512 strain realizations of the synthetic field. 

In the second numerical experiment, the strain matrix components were obtained from 
a 512 3 spectral DNS of the Navier-Stokes equations at Reynolds number R\ ~ 200. The 
strain time series were recorded along 512 fluid paths. The r.m.s. of the different strain 
components in this case ranged from 6.6 to 9.3 and the correlation time was approximately 
0.08. Thus, the correlation time in this case is of the same order as the inverse strain rate 
which is natural for the real Navier-Stokes turbulence where both values are of order of 
the eddy turnover time at the Kolmogorov scale. 

Because of the fairly short correlation time and because of the Gaussianity, the first 
numerical experiment is closer to the Gaussian white-noise analytical model than the 
second experiment where the strain is not Gaussian and has long time correlations. In a 
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Figure 1. Energy spectrum in the case of the synthetic Gaussian strain with a finite 

correlation time 

sense, comparison of results of the first experiment and the analytical solutions may be 
considered a performance test of the numerical method. On the other hand, comparison 
between the results of the first and the second experiments allows us to establish their 
sensitivity to the strain statistics. 

To compute Ij2.7|) we used a second order in time Runge-Kutta scheme with a time step 
of 20 times less than the correlation time for the synthetic case. For the simulation with 
strain from DNS, a time step identical to the DNS have been used (i.e. 200 times less than 
the correlation time of strain) In both numerical experiments, for each strain realization 
we consider a distribution of wavepackets (2048 in the synthetic case and 8192 for the 
DNS) with initial k chosen randomly in a sphere of radius |k| ~ 2. For each wavenumber 
k, the two Fourier components of the velocity U\ and u^, are chosen randomly such that 
uj = /3jje t27raj where otj, (3j are uniform random numbers in the interval [0,1] and 
7 = 2.e — 04 is a constant (function of the total number of particles) . The last component 
113 is deduced from k • u = to respect the incompressibility. However, the minimum 
value of I&3I is limited by the condition that when calculated 113 appears to be excessively 
large such wavepacket is discarded. The viscosity v was set to 10~ 6 in both numerical 
experiments. 

Figures ^and ID show the energy spectrum at several different moments of time in the 
first and the second numerical experiment correspondingly. In both cases one can see an 
excellent agreement with the theoretical -2.5 slope for time less than ta which is about 
4.5 and 0.7 for the first and the second experiments respectively. At t — td, turbulence 
reaches the dissipative scale and the spectrum accepts a log-corrected shape. 

Figures |3 and 0] show the time growth of the total energy and the energy spectrum at 
several fixed wavenumbers in the first and the second numerical experiment respectively. 
One can see that the growth is approximately exponential as predicted by the theory. 
It is interesting that the theory also predicts a change to a slower exponential growth 
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Figure 2. Energy spectrum in the case of the strain obtained from 512 3 DNS 

of the total energy when time cro sses td- Similar ef f ect is called the dissipa tive anomaly 
in the kinematic dynamo theory ijKazentsevI l|l968|) : Ichertkov et all l|l999|) ). This effect 
is consistent with figure |3] which shows that at t = td ~ 4.5 the slope for the total 
energy gets smaller and becomes approximately equal to the slope of the individual k- 
components (as predicted by the theory). For the second experiment the growth is also 
approximately exponential and the change of slope occurs at t = td » 0.7. 

Figures El and show the spectrum of the polarization at several fixed wavenumbers in 
the first and the second numerical experiment respectively. Similarly to the white-noise 
strain, in both simulations the polarization sharply decreases in time and it tends to an 
independent of k spectrum in the inertial range (with a log-correction for t > td). 

Figures d and |S1 show the spectrum of the flatness at several fixed wavenumbers in the 
first and the second numerical experiment respectively. Similarly to the white-noise pre- 
dictions, the flatness is growing with time. However, the theoretical 3/2 slope is observed 
neither for the synthetic nor for the DNS strain case. The fact that these deviations are 
observed in the same way for both the DNS strain and the synthetic strain (which is short 
correlated and therefore quite close to the white noise process) might indicate a failure 
of the numerical method to reproduce some features of higher correlators rather than a 
true deviation due to the differences in the strain statistics. Further discrepancy arises 
for the DNS strain case at very large time when the flatness is observed to decrease. 

9. Conclusion 

In this paper we introduced a description of the small-scale Navier-Stokes turbulence 
with much larger scales, the SDT model. Such nonlocal interactions dominate in 2D 
turbulence at large time. In 3D, they were shown to be responsible for intermittency 
in numerical experiments ijLaval et all l|200ll) ). The SDT model assumes that the large 
scales are Gaussian and short correlated in time and it describes turbulence in terms of 
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Figure 3. Time growth of the total small-scale energy (solid line) and of the energy spectrum at 
several fixed k (dashed lines) in the case of the synthetic Gaussian strain with a finite correlation 
time 




Figure 4. Time growth of the total small-scale energy (solid line) and of the energy spectrum 
at several fixed k (dashed lines) in the case of the strain obtained from 512 3 DNS 



the one-point /c-space correlators using the method introduced in lNazarenko et alj 12003) 
for the kinematic dynamo problem. We studied both 2D and 3D turbulence, the 2D case 
being equivalent to the problem of 2D passive scalars in the Batchelor regime. In 2D, we 
found steady state solutions for correlators of all orders. These solutions correspond to 
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Figure 5. Polarization spectrum in the case of the synthetic Gaussian strain with a finite 
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Figure 6. Polarization spectrum in the case of strain obtained from 512 3 DNS 

forced turbulence and they describe cascades of the energy and enstrophy series of invari- 
ants (two invariants at each order). The energy cascades are non-intermittent: initially 
Gaussian turbulence at the forcing scale remains Gaussian and scale invariant through- 
out the inertial range. On the contrary, the scale invariance and gaussianity break down 
for the enstrophy cascades and a regime dominated by thin elliptical structures develops 
in the fc-space. In 3D, we found that the steady state does not exist in SDT and the total 
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Figure 7. Flatness spectrum in the case of the synthetic Gaussian strain with a finite 

correlation time 




Figure 8. Flatness spectrum in the case of the strain from 512 3 DNS 



energy grows in a dynamo-like fashion. To have a realistic description of the steady state 
of the Navier-Stokes turbulence one have to compliment SDT with a model for the local 
scale interaction which will be done in future work. However, the study of the purely 
nonlocal interaction presented in this paper reveals several interesting effects which are 
likely to persist in some form when the local interactions are taken into account. In 
particular, the nonlocal interactions are shown to lead to an interesting turbulent state 
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in which all modes have plane polarization in the fc-space. Statistics of such turbulent 
fields is far from being Gaussian (in which a ll polarizations are pres ent). Similar effect 
for the dynamo magnetic fields was found in lNazarenko et al\ i|2003j) . The fc-space mo- 
ments allow also to quantify the dominant coherent structures in the fc-space which are 
responsible for intermittency in similar way as the i-space moments capture the fc-space 
structures. Note that singular structures in the fc-space are not necessarily singular in the 
x-space (e.g. a field periodic in some direction is a singular ID line in the fc-space). For 
the Navier-Stokes turbulence, the intermittent fc-space structures are found to be very 
elongated ellipsoids centered at the origin. These structures leave their signatures on the 
spectral flatness and on the scalings of the higher fc-space moments with respect to the 
order n. Further, similarly to the kinematic dynamo problem the scalings indicate pres- 
ence of the log-normal statistics of the small-scale velocity fields. Naturally, putting back 
the local interactions into the model will weaken the tendency to form elongated ellip- 
soids and plane polarized wavepackets, but some deviation from the Gaussian statistics 
in the direction predicted by these tendencies should be expected in real Navier-Stokes 
turbulence. 

We studied numerically the SDT model i|2.7[l for a synthetic Gaussian strain with 
a finite correlation time and for a strain field obtained from a 512 3 spectral DNS of 
the Navier-Stokes equation. The results show that most predictions of the white-noise 
theory, such as e.g. the energy spectrum shape, the polarization decrease and the increase 
of the spectral flatness, are also observed in these two cases of the strain. Thus, the SDT 
equations obtained for the Gaussian white noise strain are likely to be a good model for 
the nonlocal interactions in real Navier-Stokes turbulence. 



10. Appendix 

Our aim here is to derive a closed equation for the generating function Z starting with 
equation (|3.5|l . The last term in this equation is the easiest one, 

-2^fc 2 ((A|u(k)| 2 + au(k) 2 + (3u{kf)E) = -2vk 2 VZ (10.1) 

where T> is a differential operator defined in 1|3.8|) . The correlators containing factor cr^ 
can be found using the Gaussian integration by parts. In particular 

g£j 

(<TijE) = fi(- — ) = CI [\{(T m jjU m + T m ,ijU m )E) + 2a{T m ,ijU m E) + 2(3(T m ^u m E)\ , 
0O~ij 

(10.2) 

where we have used the definition l|3.ti[) . Here, ^ m ,ij is a response function, 



r m ,ij = (io.3) 

OCT,; 



Differentiating l|2.7|l with respect to ct^ and using whiteness of the strain tensor we get 
r m ,ij = {h dj + ^-{l-ki di))u m + ( 2k ™ kl - $mi)uj. (10.4) 
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In what follows we will use the turbulence isotropy, in particular, expressions of type 



{(uiUj +TijUi)E} 
(uiUjE) 
(HiUjE} 



1 ' 'u| 2 £>(% - ^) 



(d-1) 
1 

1 



Substituting I|1U.4|I into l|10.2[l and using the above isotropy relations we have 

M = WW - 5 -fhdi)z + - 5 -f)vz 

where V is a differential operator defined in (|3.8[1 . This allows us to find the first term 
on the RHS of ij!T5|) . 



(10.5) 
(10.6) 
(10.7) 

(10.8) 



kidj (aijE) = Q 



^ l ^k 2 Z kk + - 1 {2V + d z -l)kZ k + 2VZ 



1 

d [ 



(10.9) 



Similarly, the other three terms on the RHS of (|3.5|l can be obtained by the Gaussian 
integration by parts and using the response function (|10.4|) and the isotropy. After a 
lengthy but straightforward algebra one gets 



Hcr m i(u m ui+uiu m )E) = -2\tt 



and 

2a(a m iu m uiE) = 2aQ 



(d - - d )Z x + 2(1 - ^)VZ X + \(Z afJ - Z\\) + - d kAZ x 

(10.10) 



2 2 1 

(- - d)Z a + -VZ a - -kidiZ a - 2XZ\ a - 2/3Z A a - 2aZ aa 
d da 



(10T1) 

The 4th term can be obtained from l|10.11|) via interchanging a with j3 and u with u, 



2(3(a m iu m uiE) = 2pfl 



{^-d)Z p + ~VZp - ~hdiZ - 2\Z xp - 2aZ\\ - 2[3Z^ 



(10.12) 

Putting expressions lfl"tI3|) . (|10.10|l . I|10.11|) . (|10.12|l and ifTTH)) . we have the following 
final equation, 

Z = fi[(l- -)k 2 Z kk + -UV + d 2 - l)kZ k + 2(l - - + d)VZ - -V 2 Z 
1 a a ad 

+2(A 2 + 4af3)Z xx + 2X 2 Z a0 + 8XaZ aX + 8\f3Z fJX + 4a 2 Z aa + 4/3%^] 
-2uk 2 VZ, (10.13) 
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